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We extend the Krori-Barua analysis of the static, spherically symmetric, Einstein-Maxwell field 
equations and consider charged fluid sources with anisotropic stresses. The inclusion of a new 
variable (tangential pressure) allows the use of a non-linear, Chaplygin-type equation of state with 
coefficients fixed by the matching conditions at the boundary of the source. Some physical features 
are briefly discussed. 
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I. INTRODUCTION 

A major issue of the static, spherically symmetric Ein- 
stein field equations in GR is to find out interior solutions 
which are free of singularity. It has shown that an un- 
charged incompressible fluid sphere of mass m cannot be 
held in equilibrium below certain radius a = ^p and even 
demands a larger value for a related to physically reason- 
able equation of state (EOS) [l|. Regarding stability of 
the model Stettner [2] argued that a fluid sphere of uni- 
form density with a net surface charge is more stable than 
without charge. Therefore, a general mechanism have 
been adopted by the relativists to overcome this singular- 
ity due to gravitational collapsing of a static, spherically 
symmetric fluid sphere is to include charge to the neu- 
tral system. It is observed that in the presence of charge 
either gravitational attraction is counter balanced by the 
electrical repulsion in addition to the pressure gradient 
[3| or inhibits the growth of space-time curvature which 
has a great role to avoid singularities Q. According to 
Ivanov [5[ the presence of the charge function serves as 
a safety valve, which absorbs much of the fine tuning, 
necessary in the uncharged case. However, in connection 
to this we would like to mention here a special kind of 
mechanism to avert singularity as used by Trautman [6] 
under Einstein-Cartan theory where physical entity spin- 
torsion supposed to act as an agent of repulsive effect. 

A large amount of works on charged fluid spheres are 



available in the literature (an exhaustive discussion with 
various classification schemes regarding sources for the 
Reissner-Nordstrom (RN) space-time can be obtained in 
Ref. [5]). However, in connection to singularity wq would 
like to mention here that Efinger 7] , Kyle and Martin Q 
and Wilson Q have found relativistic internal solutions 
for static charged spheres, but none of these solutions 
is absolutely free from singularities. On the other hand, 
spheres of charged dust have been investigated by Bonnor 
1C|, Bonnor and Wickramasuriya [11] and Raychaudhuri 
12j . Among all these investigations it is observed that in 
Efinger's solution the metric has a singularity at the ori- 
gin (r = 0) whereas the solutions due to Kyle and Martin 
[8| and Wilson Q do not have any interior singularities. 
However, it is argued by Junevicus [13| that in both the 
above cases the metrics may have singularities at points 
other than the origin so that restrictions have to be im- 
posed on the sphere to avoid them. According to him 
the fluid sphere solutions of Kyle and Martin [8j , Wilson 
9J, Kramer and Neugebauer [14| and Krori and Barua 
15j are of special interest since, with the imposition of 
suitable conditions, they are completely free of metric 
singularities and satisfy physical considerations (fo r the 
discussion and analysis of stability vide Refs. [16| and 
17J respectively). 
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We would like to note here, specially, the works of Krori 
and Barua (KB) [15[ and Junevicus [13[ which are the ba- 
sis of our present investigation. Krori and Barua (KB) 
[15| constructed static, spherically symmetric solutions 
of the Einstein-Maxwell equations based on a particular 
choice of the metric components goo and gn in curva- 
ture coordinates. Assuming that the source is a charged 
fluid with isotropic stresses, the three independent Ein- 
stein equations were reduced to linear algebraic equations 



for the energy density p(r), pressure p(r) and the square 
of the electric field, E(r) 2 . In this approach the inde- 
pendent Maxwell equation is used to obtain the charge 
density from the pre-determined form of E(r). A special 
feature of KB [15] solutions is that they are singularity 
free. A thorough analysis of this singularity-free KB |15| 
solution has been done by Junevicus j!3| . The main as- 
pect of his investigation is to fix up the constants involved 
in the KB [15| metric in terms of the physical constants 
of mass, charge and radius of the source. He also in- 
vestigated the conditions for physical relevance leading 
to a functional dependence of the ratio of mass-to-radius 
on the ratio of charge-to-mass and also to upper and 
lower limits on these ratios. In his recent work on static 
charged perfect fluid spheres in general relativity Ivanov 
Q has also observed that the solution of KB 15] which 
are fixed by the junction conditions is non-singular and 
the positivity conditions are satisfied. 

In connection to the above discussion on KB solution 
|15j it is to be mentioned regarding the very recent work 
on these solution by Varela et al. [18[. The work deals 
with self-gravitating, charged and anisotropic fluids to 
solve the Einstein-Maxwell equations. In order to discuss 
analytical solutions they 15 1 extend KB method |T5| to 
include pressure anisotropy and linear or non-linear equa- 
tions of state. The obtained solutions satisfy the energy 
conditions of general relativity and have the following fea- 
tures: (1) spheres with vanishing net charge contain fluid 
elements with unbounded proper charge density located 
at the fluid-vacuum interface; (2) inward-directed fluid 
forces caused by pressure anisotropy may allow equilib- 
rium configurations with larger net charges and electric 
field intensities than those found in studies of charged 
isotropic fluids; (3) links of these results with charged 
strange quark stars as well as models of dark matter in- 
cluding massive charged particles are possible, and (4) 
the Van der Waals equation of state leading to matter 
densities constrained by cubic polynomial equations is 
considered. 

Our present investigation of static, spherically sym- 
metric charged fluid sphere distribution is in continua- 
tion of the above work of anisotropic fluid source [18j . 
This means that the radial and tangential pressures are, 
in general, unequal so that the simplest relation between 
them may be assumed as pt = np r , (n ^ 1) 19]. How- 
ever, there are several other forms of anisotropic rela- 
tionship between pressures can be noted in the literature, 
e.g., in connection to electromagnetic mass model Her- 
rera and Varela [20| introduced a condition of anisotropy 
in the form p t — p r = gq 2 r 2 where g is a constant 
having non-zero value whereas Barreto et al. [2l| de- 
fine the degree of local anisotropy induced by charge as 
Pt — Pr = t~ j where E is the local electric field inten- 
sity, to consider self-similar and charged radiating fluid 
spheres as anisotropic sources. 

Recently, scientists show great interest on Chaplygin 
gas EOS in order to explain accelerating phase of the 
present Universe as well as to unify the dark energy and 



dark matter. As Chaplygin gas EOS is a specific form 
of polytropic EOS so it looks promising to describe dark 
energy spherically symmetric charged obj ects, generally 
termed as dark stars in the literature [22H25J . As a pos- 
sible mechanism of formation it is argued by several in- 
vestigators that the first stars to form in the Universe, 
at redshifts z ~ 10 — 50, may be powered by dark matter 
annihilation for a significant period of time rather than 
nuclear fusion [2q - l29l |. On the other hand, it is believed 
that dark energy exerts a repulsive force on its surround- 
ing and this repulsive force, likewise electric charge, may 
prevent the star from gravitational collapse. Therefore, 
people have speculated that a massive star does not sim- 
ply collapse to form a black hole, instead to the formation 
of dark energy star with a final configuration without nei- 
ther singularities nor horizons J30H33l |. 

However, among the above mentioned dark star mod- 
els we are specially interested to the work of Bertolami 
and Paramos |22f where they, like us, have used the gen- 
eralized Chaplygin gas (GCG) EOS in special reference 
with anisotropic pressure though our motivation and ap- 
proach to solve the spherically symmetric gas model is 
quite different from them. The scheme of the present 
investigation is therefore as follows: we write down the 
four independent Einstein-Maxwell equations. By allow- 
ing the radial (p r ) and tangential {pt) pressures to be 
different, we have found out the six variables p r , pt, p, 
e = E(r) 2 , A, v where the other parameters are, respec- 
tively, matter-energy density, electric field intensity and 
metric potentials (Sec. II). Ivanov [5J has explained the 
usual difficulties that generally arise when we combine 
equations of state (even a linear one) with the field equa- 
tions. Interestingly, we are here dealing with a non-linear 
EOS and are able to find solutions using an algebraic 
method (we do not solve differential equations). Adding 
the non-linear Chaplygin gas EOS p r = Hp — — (where 
H and K are two arbitrary constants) and using the KB 
ansatz for A and v we get four algebraic equations for 
p, p r , p t , e = E(r) 2 . Using the independent Maxwell 
equation we determine the charge density a from e (Sec. 
III). We present and discuss the necessary matching of 
the solutions and the related boundary conditions in the 
Sec. IV which allow us to find out the expressions for H 
and K with their physical features through the graphi- 
cal plots. Also, energy conditions have been discussed in 
detailed (Sec. V). Some concluding remarks are made in 
Sec. VI. 



II. BASIC EQUATIONS 

The KB [15J metric is given by 



ds 2 



3^dt 2 + e x ^dr 2 +r 2 (d0 2 



'■) (1) 



where A(r) = Ar 2 and v(r) = Br 2 + C with arbitrary 
constants A, B and G. 

The most general energy momentum tensor compatible 



with spherically symmetry is 

T£ = {p + PrWUv - Vrgt + (pt - Pr)rfriv 

with 

u^ Ufl = -rfri/j, = 1 

The Einstein-Maxwell equations are 



A' 1 



1 



= 8tt P + E 2 



(2) 



(3) 



where, as mentioned earlier, A, B and C are some ar- 
bitrary constants. It is of interest to note that these 
constants were determined by Junevicus [13] in terms of 
the physical quantities mass, charge and radius of the 
source. 

Along with the above ansatz let us also use generalized 
Chaplygin gas EOS for the charged fluid as [34[ 



p 



(9) 



- — = 8ir Pr - E 2 



(4) 



where H and K are two positive constants. 
Equations © and (H() implies 



1 

r 

and 



\{v') 2 + v"-\\'v' + \(v'-\') 



Snpt + E 2 

(5) 



(r 2 E)' = 4irr 2 <je^. (6) 

Equation (6) can equivalently be expressed in the form 

q{r) 



r Jo 



<±nr ae 2 dr = 



(7) 



where q(r) is the total charge of the sphere under 
consideration. 



(p + Pr ) s/ (r) = ^±^ e -^. (10) 



-in 



From equations ([9]), ((H) and ((T0|) . we get the following 
solution set: 



P = 



f+^f 2 + 2hK 



(11) 



III. SOLUTIONS 

Now, we consider KB ansatz: 

X(r) = Ar 2 , v{r) = Br 2 + C 



(8) 



Pr = / 



7 + V/ 2 + 2fex \ 



(12) 



Pt= 8n 



-Ar 2 



2e- Ar (B - A){2 + Br 2 ) - — 1 - e~ Ar 



8tt 



/ + V/ 2 + 2hK 



(13) 



E 2 = 2Ae- Ar ' 2 + 4(1 - e ~ Ar2 ) - ^ 



f + v / f 2 + 2hK 



(14) 



q 2 =2Ar 4 e- Ar +r\l-e~ Ar ) 



/ + y/p + ZhK 



(15) 



where h — 2(1 + H). We observe that for finite values 
of the physical parameters h ^ 0, so that H ^ — 1. 



IV. BOUNDARY CONDITIONS 



The RN metric [35j, |36( is given by 



ds 2 



2m Q 2 



dt 2 



2m 

r 



dr 2 +r 2 (d6 2 + sin 2 



(16) 



r 



To match our interior metric with the above exterior one 
we impose only the continuity of g u , g r r an d -f^ 1 across 
a surface, S, at r = a. This yield the following equations: 



_ 2m Q_ _ Ba 2 +C 

1 ~t~ o — C , 



l-^ + ^=e- 

a er 



^2 



Bae 



Ba^+C 



(17) 



(18) 



(19) 



a" a J 
Therefore, from the above equations, one can easily get 



a z 



1 



B 



C = ln 



m 
a 

2m 
a 



Q 2 



Q 2 



2m Q 2 

— + -\ 

a a z 



2m Q 2 



i 2m 



Q 2 



(20) 



(21) 



(22) 



Note that no extra assumption on the value of e(r) at 
r = a is required here and we therefore obtain e(a) — 

z*t. This result was expected as a consequence of the 
matching conditions at r = a (absence of thin shell). 
On the other hand, the electric field at r — is zero. 
Therefore, the energy conditions at r = involve only 



the central value of density, as well as for the radial and 
tangential pressures. As thin shell exists, so we could not 
use boundary condition for e at a which means one can 

not get e(o) = fr. 

Let us now impose the boundary conditions 



p r (a) = 0, e(0) = E(0) = 0. 



(23) 



We obtain two independent equations which are readily 
solved for H and K as functions of the source parameters. 
We note that -E(O) = implies 



Fe 



Aa z 



(Fe Aa2 



3A 



and p r (a) = implies 

r , ~F + VF 2 + 2hK 



2hK 



(24) 



= 



(25) 



- (A±B). P -Aa 2 
Ait 



where F = f(a) 

From the above two equations (|24f and (|25j) , one could 
find the values of two unknowns H and K in terms of 
A, B and a, in other words, in terms of mass, charge 
and radius of the spherically symmetric charged objects. 
Therefore, through a simple mathematical exercise, we 
have the expressions for the constants as follows 



H 



3A(B - A) 



[9A 2 -4{A + B) 2 e- 



-2Aa 2 ] 



(26) 



K = 



(A + B)e- 

47T 



2 1 



3A{B - A) 



[9A 2 -4(A + B) 2 e~ 2Aa2 ] 



- 1 



(27) 



One can note that at r = 0, E(0) = 0, p(0) = |£, where 
Pr(0) = ipt(0) = ^§^. Also, the curve profiles (Figs 1 
- 3) for the parameters, p, p r etc. indicate no singularity 
presents inside the star. 



Fg = -Br(p + p r ). 



(31) 



V. TOV EQUATIONS 

The generalized Tolman-Oppenheimer-Volkov (TOV) 
equation as presented by Ponce de Leon 37] is 



M G {p + Pr ) 



dp r 

dr 



q a 



-iPt ~Pr) = 0, 
r 

(28) 
where Mq — Mq(t) is the effective gravitational mass 
inside a sphere of radius r and q = q(r) is given by (15). 
The effective gravitational mass is given by the expression 



Mair) = -rV 



(29) 



derived from the Tolman-Whittaker formula and the 
Einstein-Maxwell equations. 

It is important to note that the above equation de- 
scribes the equilibrium condition for charged fluid ele- 
ments subject to gravitational, hydrostatic and electric 
forces, plus another force due to anisotropy factor which 
is a measure of the pressure anisotropy of the fluid com- 
prising the charged body. Combined with (11), (12), (13) 
and (15), the above equation takes the form 



F h = 



dp r 

dr 



F P = aEe~ 



F a = - (pt -Pr) ■ 
r 



(32) 



(33) 



(34) 



F g + F h +F e + F a = 0, 



(30) 



The profiles of F g , Fh, F e and F a for sources are shown 
in Fig. 4. This figure indicates that Fh is comparatively 
small. Thus the hydrostatic force has a negligible effect 
in spite of the static equilibrium is attainable due to pres- 
sure anisotropy, gravitational and electric forces. 

Though several specific equation of states for p r (p) are 
used in literature, but a very little is known for the much 
less intuitive second equation of state pt(p)- The equation 
of state parameter u> t = — for the anisotropic object can 
be obtained directly from equations (11) and (13), which 
is given by 



El 

9 



w* = 



2e-^ 2 (B-A)(2 + i?r 2 )-^( 



1-e 



-Ar- 



8n ( /+V7W 



8tt 



( f + Vf 2 +2hK 
{ h 



(35) 



Fig. 5 shows that the variation of u>t against r. The 
equation of state for anisotropic charged fluid is positive 
and confined within < cc> t < 1 i.e. charged fluid is non 
exotic in nature. 



VI. ENERGY CONDITIONS 



It is well known for the charged fluid that the null 
energy condition(NEC), weak energy condition (WEC) 
and strong energy condition (SEC) will be satisfied if 
and only if the following inequalities hold simultaneously 
at every point within the source: 
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FIG. 1: The density parameter p is shown against r. 
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FIG. 2: The radial pressure p r is shown against r. 



E 2 


(36) 


p + p r > 0, 


(37) 


E 2 

P+Pt + —>0, 

47T 


(38) 


.B 2 

/3 + J?r + 2p t + — - >0, 

47T 


(39) 



Direct plotting of the left sides of (9)-(12) show that 
these inequalities are satisfied as well at every r (see Fig. 
6). 

At this point we feel it is required to determine whether 
specific choices of mass, charge and radius lead to solu- 
tions satisfying the above energy conditions at r = 0. 
A close observation of the equations ([2l7| and ([2"TT) sug- 
gest us to adopt here some adimensional quantities which 
can be defined as a = a 2 A and j3 = a 2 B. We re- 
strict our attention to solutions satisfying a > 0, (3 > 0, 
< 2a. These solutions satisfy the four energy con- 
ditions of general relativity, viz., Null Energy Condi- 
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FIG. 3: The transverse pressure pt is shown against r. 
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FIG. 4: Four different forces acting on fluid elements in static equilibrium is shown against r. 



tion (NEC), Weak Energy Condition (WEC), Strong En- 
ergy Condition (SEC) and Dominant Energy Condition 
(DEC). 

For a sphere of radius a, mass M and charge Q, equa- 
tions (2.13) - (2.15) in [H| (equations flU} - flU) in our 
present case) are alternatively expressed in the adimen- 
sional forms 



• In (1 - 2/i + X z 



1 - lii + x 2 



(40) 



(41) 



C = In (1 - 2/i + X 



[t-xr 



1 - 2/i + x 2 ' 



(42) 



IOI 



where /t = — ,x — ■ It is very important that the 
field equations can eventually be expressed in terms of 
these adimensional constants, the adimensional variables 
p = a 2 p,p r = a 2 p r ,pt — a 2 pt,e = a 2 e, and the adi- 
mensional radial coordinate x = -. We have seen that 

a 

particular values of the adimensional parameters p, x de- 
termine the adimensional KB constants a, (3, C which in 
turn determine A and v at every x € [0,1]. The values of 
p and x are restricted by the condition that no horizon 
is included in the external region described by the RN 
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FIG. 5: The variation of equation of state parameter uit is shown against r. 
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FIG. 6: The variation of left hand side of the expressions of energy conditions are shown against r. 



metric. 

We consider all possible roots of the equation goo = 0. 
The radius of the charged sphere a is big enough so that 
no horizons are included in the external RN metric. The 
three possible cases are: 



1 -L 2 

satisfies x < M < —? 



A. Two real roots 



B. One real root 



H 2 > x 2 - 



v = x- 



We choose 1 > /j, + y/jj, 2 — x 2 - Therefore, x < 1 an d H We choose 1 > \x = X- 



C. No real roots 

fj, < x, otherwise arbitrary. 

The selected values of p, and \ determine values of a 
and /3 which should satisfy the energy conditions. An- 
other acceptability condition is that e(x) > for every 

xe [0,1]. 

The arising expressions for p and p r and j>t can be 
evaluated at r = 0. Hence we find the energy density 
and pressures at r = as 



3A 



3a 



PrO = ^PtO 



2B- A 2/3 - a 



8n 



87TO 2 



(43) 



(44) 



Now, 


EOS at 


r = 


is 

PrO _ 
PO 


- rr 


, 








m 


2/3 -a 
3a 


< 


/3 
2a 


< 1 



(45) 



(46) 



Notice that p and p r are decreasing functions of r 
(these can be shown by plotting the graphs of p and p r 
or one can find ^ < and -J^f < i.e. p and p are 
decreasing functions of r). Since, at r = 0, they assume 
fixed values, po and p r o- So, p and p has a maximum at 
r = 0. We have checked that p = 0, pj, = and p' ' < 0, 
P" < 0. 



Now, The energy conditions |38( at the centre can be 
written as : 



VII. STABILITY 



(i) NEC: p r o + po > =*> a + /3 > 

(ii) WEC: p r0 + p >0^a + ,9>0&p >0^ 
a> 

(iii) SEC: p rO + po>0^a + /?>0& 3p r0 + po > 
^/3>0 

(iv) DEC: po >| Pro \ ^ 2a > /3 > a 

The characterization of dark energy fluid is the viola- 
tion of one of the SEC, more specifically, that one related 
the Raychaudhuri equation 25, 39] . If the second of the 
WEC is violated, we have a phantom dark energy fluid. 



Bertolami and Paramos [22| argue that if one assumes 
that the GCG tends to a smooth distribution over space 
then most density perturbations tend to be flattened 
within a time scale related to their initial size and the 
characteristic speed of sound. 

One of the important "physical acceptability condi- 
tions" for anisotropic matter are the squares of radial 
and tangential sound speeds, defined by 



dp r 
dp 



H 



K 



f+^/P+2hK 
h 



(47) 



and 



J 



dp t 
dp 



y// 2 + 2hK 
AprA(A + B) 



4rA(B- A)(l + Br 2 ) 



2A 



(48) 



should be less than the speed of light [40J, |41| 



From the above equation (J47J) , an important aspect can 
be observed that the squared of radial sound velocity is 
always positive irrespective of matter density and hence 
this is always positive even in the case of exotic matter. 
The Figs 7 and 8 show that these parameters satisfy the 



inequalities < v sr < 1 
within the charged fluid. 



and < Vg t < 1 everywhere 



Now, we use Herrera's approach [40J to identify poten- 
tially unstable or stable anisotropic matter configuration 
known as the concept of cracking (or overturning). Since, 
< vL. < 1 and < w 2 t < 1, therefore, according to 
SS E| , | v 2 st - v 2 sr |< 1. The Fig. 9 of the model also 



supports this. 
Now, 



implies 



-1 <vl-vt <1 



-Kvl-vt 



<0 , 



< v 2 st - v 2 sr < 1 



potentially stable 



potentially unstable 



One can note that the region for which v 2 t < v 2 r is po- 
tentially stable region and the region for which v 2 t > v 2 r 
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FIG. 7: The variation of radial sound speed v gr is shown against r. 
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FIG. 8: The variation of tangential sound speed vlt is shown against r. 



is potentially unstable region. If v 



keeps the same 



sign everywhere within a matter distribution, no crack- 
ing will occur. The curve profile (Fig. 10) for v* t — uf r 
indicates that there is a change of sign and thus alter- 
nating potentially unstable to stable region within the 
distribution. 



VIII. MINIMUM MASS-RADIUS RELATION 

From, the above analysis indicates that our model is 
very much unstable within radius 1.5 unit. But, the the 



configuration is stable within 1.5 < r < 8. In a recent 
paper, Andreasson [421 ] has discovered a surprising result 
as: 



M < 



R 



91 

3R 



for a lower bound on the radius R of a charged sphere 
with mass M and charge Q. 

The inequality is shown to hold for any solution which 
satisfies p r + 2p t < p. 

The plot (Fig. 11) for p r + 2p t — p against r indicates 
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FIG. 9: The variation of | v st — v ar | is shown against r. 
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FIG. 10: The variation of Vg t — v%, is shown against r. 



that in the region 1.5 < r < 8, p r + 2p t — p is nega- 
tive. Since, our model is stable within 1.5 < r < 8, so 
Andreasson's relation holds good for our model. 

It would be interesting to made some comments re- 
garding opposite situations with the maximum mass- 
radius relation. By using the static spherically symmetric 
gravitational field equations Buchdahl [l| has obtained 
an absolute constraint of the maximally allowable mass- 



radius for isotropic fluid spheres of the form ^- < | (for 
a generalized expression we refer the work of Mak et al. 

It is worthwhile to calculate effective gravitational 
mass which due to the contribution of the energy den- 
sity p of the matter and the electric energy density -^ 
and can be expressed as 
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FIG. 11: The variation of 2pt + p r — p is shown against r. 
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r 2 dr = -R 



Va 



,AR 2 )- 



2VA 



-,AR' 



(49) 



where 7 (|, AR 2 ) is the lower incomplete gamma func- 
tion. In Fig. 12, we plot the mass-radius relation. We 
also plot effe cU v= a g a i ns t R (see Fig. 13) which shows 

that the ratio — e// ^ ct '" e is decreasing even if the radius 
is increasing with the mass. 

According to Ponce de Leon [37| the energy conditions 
require 

e~ A < 1 , v' > 
. These relations lead to a maximum charge is as follows: 

*>£ 

~ M 

where M and q represent the total mass and charge of 
the charged sphere of radius R. 

Assuming R = 8 Km and -^ = 0.5, we find from equa- 
tion (15) as q 2 = 0.20372. Thus our model satisfies the 
Ponce de Leon's condition [371 ] . 

In this connection we add that the mass-radius-charge 
relation for compact astrophysical objects plays an im- 
portant role in many physical processes. The strong grav- 
itational field due to the density of the matter inside 
the stars indicates that a strong electric field due to the 
electric charge is possible to exist. The effect of electric 
charge in compact stars assuming that the charge distri- 
bution is proportional to the mass density were studied 
recently by several authors [la . l44| - |4q] . 



IX. CONCLUSIONS 

In this paper we are checking the energy conditions 
only at the center of the charged sphere. It would be 
convenient to extend the analysis to other points within 
the sphere. A scries method like the one used in equations 
(32) and (34) by KB [l| might be useful. 

Unlike the work of Bertolami and Paramos [22| 
where they have used the generalized Chaplygin gas 
(GCG) EOS in special reference to anisotropic pres- 
sure we generate the solutions for KB metric under 
Einstein-Maxwell space-time. So a natural question - 
Does the result go over into the solution of Bertolami 
ct al. for isotropic stresses? The straight forward 
answer is No. This is because we have extended 
KB approach assuming singularity free form of the 
metric ansatz to charged anisotropic source with non 
linear, Chaplygin type equation of state. Therefore, 
whether our solution corresponds to Chaplygin dark 
star needs a special verification, specifically whether 
charged Chaplygin dark star does exist demands further 
investigation. In a similar fashion one may raise the 
question that does the result go over into an exact 
solution for an EOS p = Hp (K = 0)? The answer 
this time is affirmative, as our solution coincides with 
the solution obtained by Varela et al. [18| for an EOS 
p = Hp (K = 0). One can see easily that our results 
go over into the expressions obtained by Varela et al. 
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FIG. 12: The variation of M e ff ec tive is shown against R. 
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FIG. 13: The variation of 



is shown against R. 



with a\ — 0. Also, It may be interesting to extrapolate 
the present investigation to the astrophysical bodies, 
specially quark or strange stars with radius around 8 km. 



APPENDIX: ANALYSIS OF JUNCTION 
CONDITIONS 

We can note that the metric coefficients continuous at 
the junction i.e. at S where r = a. However, this does 
not mean that the metric coefficients be differentiable at 



the junction. The affine connection may be discontinuous 
there. The above statement can be quantified in terms 
of second fundamental form of the boundary. 

The second fundamental forms associated with the two 
sides of the shell are [47J45CJI 



K ?:=-< 



d 2 x v 



v dX a dX$ 



(A.50) 



where njj are the unit normals to S and can be given by 



nt =± 



a/3 



df df 



dX a dXP 



df 

dX v 



(A.51) 
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with n M n M = 1. Here, £* are the intrinsic coordinates 
on the shell with / = is the parametric equation of 
the shell S and — and + corresponds to interior (our) 
and exterior (RN). It is to be noted that since the shell 
is infinitcsimally thin in the radial direction there is no 
radial pressure. Using Lanczos equations J47H50l |. one 
can find the surface energy term S and surface tangential 
pressures pg = pj, = pt as 



1 

Ana 



(A.52) 



Pt 



8ttcl 



av 
~~2 



-i + 



(A.53) 



The metric functions are continuous on S, then one 
finds 



2 = 0, 



(A.54) 



Pi 



1 
8ira 



M 
a 



2M 

a 
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a 2 



Aa\ 1 



2M 

a 
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a 2 



(A.55) 



r 



Hence one can match our interior solution with an 
exterior RN solution in the presence of a thin shell. The 
whole space-time is given by our metric and RN metric 
which are joined smoothly. 
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